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Abstract 

The  CO2  poisoning  effect  on  anodes  for  the  proton-exchange-membrane  fuel  cell  (PEMFC)  was  examined  by  model  studies.  It  was 
assumed  that  the  reverse  water  gas  shift  reaction  (RWGS)  is  the  origin  of  the  CO2  poisoning  effects.  The  relation  between  the  anode 
polarisation  losses  and  the  catalytic  properties  of  the  catalyst  was  investigated  with  a  kinetic  model  and  with  a  fuel  cell  model  including 
finite  utilisation  of  the  fuel.  It  was  found  that  the  main  effect  of  the  occurrence  of  the  reverse  water  gas  shift  reaction  is  that  a  large  part  of 
the  catalytic  surface  area  becomes  inactive  for  hydrogen  dissociation.  Desorption  of  CO  formed  by  reduction  of  CO2  followed  by  transport 
down  the  anode  gas  channel  and  subsequent  re-adsorption  on  the  catalyst  was  shown  to  play  a  minor  role.  In  reformate  gas,  where  besides 
CO2  traces  of  CO  are  present,  CO2  poisoning  will  have  the  largest  effect  when  the  CO  content  is  small  and  at  relatively  low  current  density. 
©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

For  a  low  temperature  proton-exchange-membrane  fuel 
cell  (PEMFC)  pure  hydrogen  is  the  preferred  fuel.  How¬ 
ever,  Ht  has  a  low  energy  density  and  at  present  there  is 
no  infrastructure  for  hydrogen  distribution  and  storage  on  a 
large  scale.  Conventional  fuels  like  gasoline,  natural  gas  or 
methanol  do  not  have  these  disadvantages.  In  order  to  be 
used  in  a  fuel  cell  they  have  to  undergo  a  reforming  process. 
In  the  reforming  process  the  fuel  is  converted  by  water  or 
air  to  a  mixture  of  hydrogen,  carbon  dioxide,  nitrogen  and 
water,  and  trace  amounts  of  other  compounds  among  which 
carbon  monoxide  is  most  important.  The  efficiency  of  the 
system  is  affected  by  this  reforming  process,  not  only  be¬ 
cause  of  the  reforming  reaction  itself  but  also  because  the 
operation  of  the  PEMFC  on  reformate  gas  is  less  efficient 
than  the  operation  on  pure  Ht. 

The  detrimental  effect  of  small  amounts  (ppm  level)  of 
CO  on  the  cell  performance  is  well  known.  The  presence 
of  inert,  non-condensing  gases  like  nitrogen  results  in  ad¬ 
ditional  dilution  effects  that  increase  with  utilisation  [1,2]. 
Carbon  dioxide,  which  is  present  in  the  order  of  25%,  also 
contributes  to  this  dilution  and  utilisation  loss,  but  it  has 
often  been  recognised  [3-5]  that  the  effects  are  larger  than 
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those  of  an  equivalent  amount  of  nitrogen,  i.e.  that  in  a  fuel 
cell  CO2  is  not  inert. 

Reduction  of  CO2  on  Pt  by  Pt  hydrides  in  acidic  solutions 
in  the  0-250  mV  versus  NHE  region  was  already  reported 
by  Giner  [6] .  The  product  was  found  to  be  a  species  resem¬ 
bling  CO  adsorbed  on  Pt.  Further  work  has  shown  that  the 
reduction  of  CO2  needs  the  presence  of  hydrogen  atoms  ad¬ 
sorbed  on  the  catalyst  surface  [4,5]  but  that  the  reaction  can 
also  take  place  as  a  normal  surface  reaction,  i.e.  not  as  an 
electrochemical  reaction  [3].  It  has  been  assumed  that  the 
reaction  taking  place  is  the  so-called  reverse  water  gas  shift 
reaction  (RWGS)  in  the  form: 

C02  +2M  —  H=^M  —  CO  +  H20  +  M 

where  M  represents  a  catalyst  site.  The  adsorbed  hydrogen 
can  be  formed  by  electrochemical  reduction: 

H+  +  e  +  M=jM-H 

or  by  a  chemical  (Tafel)  reaction: 

H2  +  2M  =)>  2M  —  H 

In  a  fuel  cell  operating  at  slightly  positive  anode  poten¬ 
tials  (order  50  mV)  and  with  hydrogen-rich  gas  supplied  at 
the  anode  side  the  chemical  pathway  will  probably  prevail. 
In  the  form  written  above  the  product  species  besides  water 
is  adsorbed  CO.  Whereas  there  is  no  doubt  that  the  prod¬ 
uct  is  a  species  adsorbing  on  the  catalyst  surface,  there  is 
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Nomenclature 

bc 

Tafel  slope  of  the  CO  electro-oxidation 
reaction  (V) 

bh 

Tafel  slope  of  the  H  electro-oxidation 
reaction  (V) 

f 

molar  flow  rate  of  the  fuel  times  2 F 
(A  cm-2) 

F 

Faraday  constant  (96,485  C  mol-1) 

jav 

average  total  current  density  in  the  cell 
(A  cm-2) 

jCO 

local  generated  current  density  from 
CO-electrochemical  oxidation  (A  cm-2) 

/H 

local  generated  current  density  from 
H-electrochemical  oxidation  (A  cm-2) 

^loc 

Ku  K2,Kt,, 

local  current  density  at  position  w  in 
the  cell  (A  cm-2) 

k4,  k5 

parameters  for  Vst 

ka 

H2  adsorption  rate  constant  times 

2 F  (A  cm-2  bar-1) 

^ac 

CO  adsorption  rate  constant  times 

IF  (A  cm-2  bar-1) 

kd 

H2  desorption  rate  constant  times 

2F  (A  cm-2) 

^dc 

CO  adsorption  rate  constant  times 

IF  (A  cm-2) 

^ec 

CO  pre-exponential  electrochemical 
oxidation  rate  (A  cm-2) 

^eh 

H  pre-exponential  electrochemical 
oxidation  rate  (A  cm-2) 

kis 

rate  constant  of  the  reverse  water  gas 
shift  reaction  (A  cm-2) 

«COi 

molar  CO2  flow  rate  times  2 F  (A  cm-2) 

PX 

partial  pressure  of  gas  X  (bar) 

R 

gas  constant  (8.31  J  K-1  mol-1) 

T 

cell  temperature  (K) 

u,  u' 

total  fuel  utilisation  and  utilisation  at 
position  w 

vst 

cell  voltage  excluding  anode  losses  (V) 

W 

dimensionless  position  along  gas 
channel  (0  =  inlet,  1  =  outlet) 

Xco 

molar  fraction  of  CO  in  the  fuel 

XH 

molar  fraction  of  H2  in  the  fuel 

xco2 

molar  fraction  of  CO2  in  the  fuel 

Greek  letters 

l 1 

anode  overpotential:  potential  difference 
between  the  electron  conducting  phase 
and  the  proton  conducting  phase  in  the 
anode  (V) 

Oh 

fraction  of  the  catalytic  surface  area 
covered  by  H 

0co 

fraction  of  the  catalytic  surface  area 
covered  by  CO 

p  molar  area  density  of  catalysts  sites 

(mol  cm-2) 


some  discussion  regarding  the  true  nature  of  the  adsorbate 
[7-9].  However,  it  seems  it  is  closely  related  to  adsorbates 
occurring  in  the  adsorption  and  electro-oxidation  of  CO  and 
methanol. 

CO2  is  not  as  strong  a  poison  for  the  fuel  cell  as  CO. 
With  Pt  electrodes  and  25%  CO2  in  H2,  Wilson  obtained  cell 
performance  losses  much  smaller  than  found  with  20  ppm 
CO  in  H2  [3].  Although  this  shows  that  CO2  poisoning  is 
a  minor  source  of  performance  losses  compared  to  CO,  it 
can  still  be  significant  in  those  cases  where  due  to  good 
fuel  processing  systems  the  CO  content  is  small  (<10ppm). 
Moreover,  the  sensitivity  towards  CO2  is  not  the  same  for  all 
catalysts.  In  our  own  laboratory  we  observed  that  on  PtRu 
the  effects  were  smaller  than  on  Pt  [5],  whereas  on  PtMo  the 
effects  were  substantially  larger  [10].  Similar  trends  were 
observed  elsewhere  [11,12].  Notice  also  that  effects  of  CO2 
cannot  be  eliminated  with  an  air-bleed  [1]. 

The  RWGS  reaction  described  above  can  in  principle  pro¬ 
vide  two  sources  of  anode  poisoning  in  a  fuel  cell.  First 
there  is  the  reduced  CO2,  that,  in  the  form  of  CO  or  a  related 
species,  blocks  the  catalyst  surface  for  hydrogen  dissociation 
and  oxidation.  In  addition,  CO  formed  by  CO2  reduction  can 
desorb  from  the  surface,  enter  the  gas  flow  and  re-adsorb  fur¬ 
ther  down-stream.  Thermodynamic  calculations  have  shown 
that,  depending  on  the  water  content  and  temperature,  CO 
concentrations  of  up  to  200  ppm  can  be  achieved  in  a  re¬ 
formate  mixture  [5,13].  In  practice  concentrations  would  be 
smaller  and  depend  on  the  gas  flow.  If  the  second  mechanism 
were  dominant  it  would  be  expected  that  a  more  CO  tolerant 
catalyst  would  also  be  more  CO2  tolerant.  The  example  of 
PtMo,  already  referred  to  above,  shows  that  this  is  not  the 
case  [10-12].  On  PtMo  anode  losses  induced  by  CO2  were 
larger  than  would  be  caused  by  the  maximum  concentration 
of  CO  that  could  be  formed  thermodynamically  [10]. 

In  order  to  investigate  further  the  mechanism  of  CO2  poi¬ 
soning,  as  well  as  the  interplay  with  the  CO  tolerance  mech¬ 
anisms  displayed  by  various  catalysts,  model  studies  were 
carried  out  for  the  anode  of  a  PEM  fuel  cell  fed  on  refor¬ 
mate.  To  this  end  we  have  combined  the  model  of  Springer 
[2],  which  describes  the  kinetics  of  H2  electro-oxidation  in 
the  presence  of  CO,  as  well  as  dilution  and  CO-oxidation 
effects,  with  the  model  for  CO2  reduction  suggested  by  Bel¬ 
lows  et  al.  [4]. 


2.  Theory 

2. 1.  Kinetics 

At  an  anode  operating  under  reformate  feed  a  couple  of 
reactions  can  take  place.  The  most  important  one  is  the  dis- 
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sociative  adsorption  and  electrochemical  oxidation  of  hy¬ 
drogen.  In  addition  CO  adsorbs  on  the  catalyst  surface  area, 
where  it  can  be  electro-oxidised  by  water.  Finally  the  CO2 
present  can  react  with  adsorbed  hydrogen.  The  model  used 
here  was  based  on  the  following  set  of  reactions: 


Reaction 

Forward 

rate 

constant 

Backward 

rate 

constant 

H2  +  2M  2M  —  H 

ka 

k(\ 

1 

2M  —  H  =►  2H+  +  2e  +  2M 

^eh 

11 

CO  +  M  CO  -  M 

^ac 

^dc 

hi 

M  -  CO  +  H20  =>• 

CO2  +  2H+  +  2e  +  M 

^ec 

IV 

C02  +  2M  —  H 

M  -  CO  +  H20  +  M 

km 

V 

Reactions  (I)  and  (II)  represent  the  Tafel-Volmer  mech¬ 
anism  for  hydrogen  oxidation.  The  CO  adsorption  and  des¬ 
orption  was  modelled  by  a  simple  Langmuir  process,  based 
upon  the  adsorption  of  one  CO  per  catalyst  site.  For  sim¬ 
plicity  the  rate  constants  of  the  reactions  (I)  and  (III)  were 
not  taken  to  be  dependent  on  the  fractional  CO  coverage  of 
the  site,  unlike  in  Springer’s  model  [2].  The  CO  oxidation 
by  water  is  represented  here  by  as  a  single  step,  although 
in  general  it  will  be  at  least  a  two  step  process  involving 
the  formation  of  adsorbed  OH  that  subsequently  oxidises 
the  adsorbed  CO.  The  CO2  reduction  reaction  is  treated  as 
suggested  by  Bellows  et  al.  [4],  i.e.  adsorbed  CO  represents 
the  CO2  reduction  product.  Notice  that  the  electrochemical 
oxidation  of  CO  (IV)  can  be  considered  as  the  backward 
reaction  of  reaction  (V). 

Like  in  Springers  work  the  rate  constants  of  the  adsorption 
reactions,  ka,  kac ,  krs  will  be  expressed  in  A/cm2  (geomet- 
ric).bar.  The  other  rate  constants  are  in  A/cm2.  The  reaction 
constants  for  the  reactions  (I),  (III),  (V)  include  a  factor  twice 
the  Faraday  constant,  as  1  mol  of  gas  is  in  principle  equiva¬ 
lent  to  two  reacting  electrons.  The  rates  of  electrochemical 
reactions  depend  on  anode  overpotential.  For  the  hydrogen 
oxidation  a  hyperbolic  sine  function  was  used,  allowing  zero 
current  density  at  zero  overpotential.  The  CO  oxidation  rate 
was  expressed  by  a  normal  Tafel  behaviour.  Throughout  this 
work  we  assume  water  saturation  of  the  gas  phase,  so  the 
water  partial  pressure  is  constant  and  can  be  incorporated 
in  kec-  When  p  represents  the  molar  area  density  of  catalyst 
sites  and  F  the  Faraday  constant  the  following  equations  re¬ 
sult  for  the  fractional  hydrogen  and  CO  coverage: 

=  k*PR^1  -  On  -  Ocof  -  kd02H  -  kTsPC 02^ 

—  2keh#Hsinh-^-  (1) 

bh 

dOco 

2  Fp— —  =  kac/?CO(l  —  On  —  #co)  —  kdcOco 
at 

+  krspCO20H  -  2kec0coexp-^-  (2) 

br 


The  current  densities  are: 

ij  0 

/H  =  2fceh0Hsmh  — ,  /co  =  2kec@coexp  — ,  and 
bh  bc 

I  —  In  +  /co  (3) 


In  the  stationary  state  the  left-hand  side  of  the  Eqs.  (1) 
and  (2)  becomes  zero  and  they  can  be  written  as: 

0  =  k„pH2(l  -On-  Oco)2  -  kdOn  -  krsPCO202 

-  2keh0nsinh5-  (4) 

bh 


kac/?CO(l  —  On) 

uco  —  - 

kacpCO  +  kdc  +  kec  exp  (i)/bc) 

+  - krsPC02 - 02  (5) 

kacpCO  +  kdc  +  kec  exp  (rj/bc) 

By  substituting  Eq.  (5)  into  Eq.  (4)  On  can  be  solved  as 
one  of  the  real  roots  of  a  fourth-order  polynomial,  with  the 
conditions: 

0  <  On  <  1 ,  0  <  Oco  —  1  >  and  On  +  Oqq  <  1  (6) 


For  further  reference  we  consider  here  two  special  cases: 


(1)  No  CO2  in  the  fuel.  This  is  the  same  case  as  treated  by 
Springer.  Eq.  (5)  reduces  to: 


Oco 


kacpCO(l  -  On) 
kac  pCO  +  kdc  +  kec  exp(??/Z>c) 


(7) 


Notice,  that  if  there  is  no  desorption  or  oxidation  of 
CO  the  final  result  of  solving  (4)  and  (7)  would  be  On 
=  0,  Oqo  =  L  If  there  is  only  desorption  of  CO  (and  no 
electrochemical  oxidation)  there  is  a  limiting  value  for 
6>co-  attained  when  0\\  =  0: 


co 


pCO 

pCO  +  kdc/k 

ac 


(8) 


and  a  corresponding  limiting  current  density: 


if  =  kaPU2 


kdc  /  kat 


kdc/  kac  +  pCO 


(9) 


(2)  No  CO  in  the  fuel.  In  this  case,  Eq.  (5)  reduces  to: 


Oco  = 


krspCXL 

kdc  +  kecexp0i/bc) 


(10) 


Also  in  this  case  without  desorption  and  oxidation  of 
CO  the  stationary  state  solutions  would  be  On  =  0, 
Oco  —  1-  Again,  if  there  is  desorption  of  CO,  and  no 
electro-oxidation  a  limiting  current  density  would  arise 
at  On  —  0.  However,  since  in  that  case  also  Oco  —  0 
(Eq.  (10))  this  limiting  current  density  is  equal  to  one 
obtained  without  CO  or  CO2: 

/h”  =  kapH2  (11) 
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2.2.  Utilisation  of  the  fuel 

There  are  several  ways  in  which  CO2  will  affect  the  per¬ 
formance  if  it  is  present  in  the  anode  feed  of  a  fuel  cell. 
Apart  from  surface  blocking,  as  indicated  above  in  the  ki¬ 
netic  model,  the  dilution  and  utilisation  of  hydrogen  lead 
to  a  reduced  partial  hydrogen  pressure  and  hence  increased 
activation  and  Nernst  losses.  In  this  respect,  the  addition 
of  CO2  is  equivalent  to  the  addition  of  N2.  Moreover,  CO 
formed  in  the  CO2  reduction  process  can  desorb  into  the  gas 
flow.  The  concentrations  that  can  be  attained  will  increase 
with  decreasing  gas  flow,  i.e.  with  increasing  utilisation  of 
the  flow.  This  CO  can  cause  poisoning  further  down  the 
channel.  In  the  present  work,  the  kinetic  model  described 
above  is  combined  with  an  integration  of  the  current  density 
over  the  flow  channel,  as  described  in  [2],  but  with  a  slight 
modification  to  accommodate  the  effects  of  the  water  gas 
shift  reaction.  When  the  partial  pressures  in  the  gas  channel 
are  a  function  of  the  position  in  the  channel,  the  current  den¬ 
sity  becomes  a  function  of  the  position  in  the  channel.  Since 
the  cell  potential  is  constant  along  the  channel,  this  means 
the  anode  overpotential  is  also  a  function  of  the  position  in 
the  channel.  The  anode  overpotential  is  calculated  by: 

RT 

riloc  =  Vceii  -  yst(/l0C)  +  T!Fln  Ph!>0C  (12> 

Zr 

The  last  term  in  this  equation  represents  the  change  in 
Nernst  potential  along  the  channel.  In  this  equation  Vst  rep¬ 
resents  the  ohmic  losses  and  the  losses  at  the  cathode.  For 
Wt  the  following  empirical  expression  is  used,  suggested  by 
Kim  et  al.  [14]: 

yst(/loc)  =  K\  —  K2ln(/loc)  -  K3Iloc  -  K4exp(-K5/loc) 

(13) 

In  this  equation  K\  represents  the  reversible  open  cell  po¬ 
tential  and  the  cathode  activation  losses  at  7loc  =  1  A/cm2, 
K3  corresponds  to  the  Tafel  slope  of  the  oxygen  reduction 
reaction  at  the  cathode  and  K3  is  a  measure  for  the  ohmic 
resistance  of  the  cell.  The  last  term  in  Eq.  (13)  with  parame¬ 
ters  K4  and  K5  represents  mass  transfer  losses  at  the  cathode 
occurring  at  high  current  densities. 

For  a  given  set  of  gas  conditions  it  is  possible  to  obtain, 
in  an  iterative  procedure,  a  set  of  consistent  values  of  rj,  In, 
Icq,  and  «co2-  The  latter  variable  represents  the  amount  of 
CO2  reacting  in  reaction  (V)  and  is  calculated  from: 

«C02  =  £rspCO20H  (14) 

Notice,  that  nco  is  expressed  in  A/cm2  (geometric)  like 
In  and  Iqo,  he.  a  factor  of  twice  the  Faraday  constant  (2F)  is 
included.  The  integration  of  the  model  over  the  channel  pro¬ 
ceeds  as  follows.  The  total  flow  is  reduced  along  the  chan¬ 
nel  as  a  result  of  the  hydrogen  consumption.  The  nitrogen 
flow  remains  unchanged,  as  does  the  sum  of  the  CO  and 
CO2  flows,  since  one  molecule  CO  reacts  to  one  molecule 
of  CO2  and  vice  versa.  The  total  water  partial  pressure  is 


kept  constant,  at  saturation  level.  The  dry  molar  flow  rate 
(i.e.  excluding  water)  at  any  point  in  the  channel  is  given  by: 

1  —  rin 

/=/in- - -  (15) 

1  -  xn 


For  the  derivatives  of  the  total  molar  flow  rate  and  specific 
molar  flow  rates  we  obtain  with  respect  to  the  dimensionless 
position  in  the  channel  (0  =inlet,  1  =outlet): 


06) 


dxnj  _  _  d.rco  / 

dw  dw 


-I  co  +  «co2 


(17) 


dv'co2  / 

dw 


=  +Ico  -  «co2 


d-YN2  /  _Q 

dw 


08) 


In  the  above  equations  the  molar  flow  rate  incorporates 
a  factor  2 F.  For  simplicity  it  is  assumed  that  the  change  of 
the  total  CO2  flow  is  small  compared  to  the  total  flow  and 
may  be  neglected.  To  calculate  the  fraction  of  CO2  in  the 
flow  the  sum  of  the  CO2  and  N2  flows  is  calculated  and  the 
fraction  of  CO2  and  N2  is  calculated  by  keeping  the  ratio 
of  the  respective  fractions  equal  to  the  inlet  ratio.  There  are 
now  only  three  variables  to  be  integrated: 


d f_  _  _j  dfii  _  ~ 7h(1  -  *h) 
dw  H  dw  / 

dyco  _  -/co  +  »co2  +  xco  1 1 1 

dw  / 

In  a  fuel  cell  the  fuel  is  usually  supplied  at  a  constant 
utilisation  or  stoichiometry,  i.e.  the  inlet  flow  rate/"  depends 
on  the  average  current  density  and  is  therefore  not  known 
until  the  average  current  density  is  known.  As  shown  in 
[2]  it  is  possible  to  perform  the  integration  over  the  local 
utilisation  u'  instead  of  over  the  channel  length.  The  relation 
between  the  two  variables  is: 


r 

Jo 


(*H  +  4no  )/"V(w)  =  /  (In  +  /co)  dw  or 
(*§  +  *&)  )/n 


dw 

du' 


/co  +  In 


(20) 


Using  Eq.  (20)  we  now  obtain  equations  for  the  derivatives 
of  the  xh  and  xco  that  can  be  integrated: 


dxH 

rin  1  vin 
aH  "r  ACO 

(1  -  xh)2/h 

du' 

1  —  vin 

1  XH 

/h  +  /co 

dxco 

vin  1  vin 
aH  aCC 

1  (1  -  xn)(Ico  -  «co2  -  *co/h) 

du' 

1  —  rin 

1  XH 

In  +  /co 

(21) 


(22) 


The  third  variable  to  be  integrated  is: 

dw  /  1  \  _  1 

du'  y  (xg  +  x“0)  fm  J  /co  +  /h 


(23) 
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Using  a  fourth-order  Runge-Kutta  scheme  xh  and  xco  are 
calculated  as  function  of  the  utilisation.  After  the  integration 
fn  can  be  obtained  from: 

/av  =  /m(x'Hn+x£0)M  (24) 

Using  this  value  of /n,  the  variable  w  can  now  be  calcu¬ 
lated  as  a  function  of  u'  from  the  integration  of  Eq.  (23). 


3.  Results 

3.1.  No  CO  in  the  fuel,  zero  utilisation 

In  order  to  investigate  the  effect  of  the  reverse  water  gas 
shift  reaction  on  the  anode  polarisation  curve,  we  first  con¬ 
sider  the  case  pCO  —  0.  In  a  fuel  cell  this  would  imply  that 
the  anode  gas  does  not  contain  traces  of  CO,  but  also  that 
the  utilisation  of  the  fuel  is  zero.  At  zero  utilisation  the  flow 
is  infinitely  high  and  therefore  the  composition  is  constant, 
i.e.  CO  formed  in  the  reverse  shift  reaction  and  desorbing 
according  to  reaction  (III)  does  not  affect  the  gas  composi¬ 
tion,  nor  does  the  consumption  of  ID  and  CO2.  The  relation 
between  Oqo  and  0]\  is  now  given  by  Eq.  (10).  Notice,  that  if 
there  is  no  oxidation  or  desorption  of  CO  (kcc  and  k&c  =  0) 
the  stationary  state  will  be  that  all  catalytic  sites  are  covered 
by  CO  ( Oco  =  1,  Oh  =  0). 

In  the  first  set  of  calculations  we  consider  the  case  where 
electrochemical  oxidation  of  CO  does  not  occur,  correspond¬ 
ing  to  the  case  of  a  Pt  electrode  at  low  (<0.2  V)  overpoten¬ 
tial.  Desorption  of  CO,  however,  was  included.  The  values 
for  the  parameters  describing  the  hydrogen  oxidation  pro¬ 
cess  (reactions  I— II)  have  been  chosen  such  that  at  normal 
fuel  cell  operating  conditions  the  anode  losses  under  pure 
hydrogen  are  minimal,  i.e.  the  limiting  current  density  is 
very  high,  and  result  in  a  realistic  hydrogen  coverage  at  low 
anode  overpotential.  These  values  are  determined  by  the  ra¬ 
tio  kd/ka.  The  parameter  ke h  was  high  enough,  so  that  the 
electrochemical  step  did  not  contribute  to  the  losses.  The 
value  for  the  rate  constant  of  the  CO  adsorption  step  was 
chosen  arbitrarily,  but  the  ratio  of  the  rate  constants  of  the 
adsorption  and  desorption  of  CO  was  chosen  such  that  a 
high  CO  coverage  results  already  at  low  CO  content.  These 
parameter  values  have  all  been  listed  in  Table  1.  In  this  set  of 
calculations  the  rate  constant  of  the  reverse  water  gas  shift 


Table  1 

Values  for  the  parameters  used  in  the  kinetic  model 


ka  (A/cm2  bar) 

40 

kd 

0.5  bar  x  ka 

&eh  (A/cm2) 

4 

bh  (V) 

0.032 

kac  (A/cm2  bar) 

20 

^dc 

3  x  10_6bar  x  kac 

ke c  (A/cm2) 

0,  3  x  10~6  or  3  x  10-5  (varied) 

be  (V) 

0.06 

krs  (A/cm2  bar) 

0,  0.002  or  0.02  (varied) 

reaction,  krs,  was  varied.  Two-non-zero  values  were  chosen, 
one  that  resulted  in  medium  coverage  of  the  catalyst  surface 
area  and  one  that  resulted  in  high  coverage  of  the  catalytic 
surface  area  (see  below). 

Eq.  (4)  with  Eq.  (10)  was  solved  for  a  range  of  values  of 
the  anode  overpotential  r],  for  a  mixture  of  PE  and  CO2.  This 
resulted  in  the  polarisation  curves  plotted  in  Fig.  la.  For  each 
of  the  current  density  values  obtained  0\\  and  Bqo  are  now 
also  known,  they  are  shown  in  Fig.  lb  and  c.  These  figures 
show  that  if  the  anode  is  operated  under  hydrogen  with  inert 
CO2  {krs  =  0)  the  hydrogen  coverage  only  decreases  slightly 
with  the  current  density,  at  least  as  long  as  the  current  density 
is  far  from  the  limiting  current  density.  When  the  reverse 
shift  reaction  can  occur  ( kTS  f  0)  the  CO  coverage  becomes 
non-zero  at  stationary  conditions,  and  0]\  drops.  The  CO 


Fig.  1.  Results  obtained  with  the  kinetic  model,  /7H2  =0.8  bar,  pCC>2 
=  0.2  bar,  pCO  =  0,  for  various  values  of  the  rate  constant  of  the  RWGS 
reaction,  krs  =  0  (A),  0.002  (I)  or  0.02  (O)  A/cm2  bar.  The  dashed  curves 
were  calculated  using  krs  =  0.02  A/cm2  bar  and  kec  =  3  x  10-5  A/cm2. 
All  other  curves  were  obtained  with  kec  =  0.  Other  parameter  values  can 
be  found  in  Table  1.  (a)  Shows  the  overpotential  as  function  of  total 
current  density,  (b)  the  fractional  coverage  of  catalytic  sites  by  hydrogen, 
and  (c)  the  fractional  coverage  by  CO. 
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coverage  dco  increases  with  increasing  rate  of  reaction  V ,  i.e. 
with  increasing  krs  as  illustrated  in  the  Fig.  lc  or  increasing 
pC02  (see  Eq.  (10)),  or  with  decreasing  k dC  (see  Eq.  (10)). 
The  lower  values  of  0]\  imply  that  a  higher  overpotential 
is  needed  for  the  same  current  density,  i.e.  increased  anode 
losses.  At  higher  current  densities  the  values  of  6qo  show  a 
small  decrease,  consistent  with  the  decrease  of  0h- 

The  first  result  of  the  water  gas  shift  reaction  therefore 
is  that  part  of  the  catalytic  surface  area  becomes  ineffec¬ 
tive  for  the  hydrogen  oxidation,  and  the  overpotential  losses 
increase.  The  amount  of  surface  area  blocked  by  CO  in¬ 
creases  with  increasing  pC02.  The  fraction  blocked  further 
increases  with  increasing  krs,  or  with  decreasing  k^c.  No¬ 
tice  that  this  finite  CO-coverage  does  not  require  a  non-zero 
value  of  pCO. 

Eq.  (10)  shows  that  according  to  the  model  the  effect  of 
CO  electrochemical  oxidation  is  similar  to  CO  desorption. 
The  only  difference  is  that  the  oxidation  rate  increases  with 
the  overpotential  and  therefore  it  mitigates  the  effects  of 
the  RWGS  reaction  at  larger  current  densities,  as  is  demon¬ 
strated  by  the  results  in  Fig.  1. 

3.2.  CO  in  the  fuel,  zero  utilisation 

When  there  is  no  CO2  in  the  fuel,  or  when  the  rate  con¬ 
stant  kIS  is  zero,  this  case  corresponds  to  the  classical  CO 
poisoning  case.  When  electrochemical  oxidation  of  CO  is 
excluded  ( kec  =  0)  then  according  to  Eqs.  (8)  and  (9)  a  max¬ 
imum  coverage  of  CO  and  a  corresponding  limiting  current 
density  occur.  In  the  case  of  10,  20  and  50  ppm  in  a  dry  gas 
with  pressure  1  bar,  the  values  for  the  maximum  CO  cover¬ 
age  calculated  with  the  parameter  values  from  Table  1  are 
0.77,  0.87  and  0.94,  respectively.  With  pU2  —  0.8  bar  the 
corresponding  limiting  current  densities  are  1.7,  0.54  and 
0.10  A/cm2.  The  curves  presented  in  Fig.  2  reflect  these  re¬ 
sults.  It  is  also  shown  that  well  below  the  limiting  current 
density  the  CO  coverage  is  smaller  than  the  maximum  value 
but  it  increases  with  the  current  density  when  0\\  shows  a 
decline  towards  zero. 

When  RWGS  reaction  is  included,  the  coverage  by  CO  is 
already  at  low  current  densities  close  to  the  maximum  value. 
This  implies  a  further  reduction  of  0]  \  and  an  increased  value 
of  the  overpotential  at  a  given  current  density.  Notice  that 
due  to  the  reduction  of  0\\  the  rate  of  the  CO2  reduction 
will  also  be  lower  than  in  the  case  without  CO.  Close  to  the 
limiting  current  density  the  CO  coverage  attains  the  same 
values  as  in  the  case  of  direct  CO  adsorption.  At  this  current 
density  (?h  — »■  0  and  the  additional  effect  of  CO2  disappears. 

In  practice,  however,  CO  will  be  oxidised  when  the  over¬ 
potential  becomes  sufficiently  large.  This  reduces  the  over¬ 
potential  both  in  the  case  with  and  without  RWGS  reaction, 
as  becomes  clear  from  comparison  of  Figs.  2  and  3.  Fig. 
3  also  shows  that  at  low  current  densities  the  effect  of  the 
RWGS  reaction  is  larger  than  at  high  current  densities.  At 
low  current  densities  the  RWGS  reaction  results  again  in  a 
CO  coverage  close  to  the  maximum  value  given  by  Eq.  (9). 


Fig.  2.  Curves  showing  the  effect  of  the  reverse  shift  reaction  calculated 
according  to  the  kinetic  model.  The  gas  composition  was  p FF  =  0.8  bar, 
pC O2  =  0.2 bar,  and  pCO  lOpbar  (A),  20pbar  (H)  or  50 p, bar  (♦). 
The  dashed  curves  with  open  symbols  were  obtained  with  krs  =  0,  the 
drawn  curves  with  closed  symbols  with  krs  =  0.02  A/cm2.  All  curves  were 
calculated  with  kec  =0.  Other  parameter  values  can  be  found  in  Table  1 . 
(a)  Shows  the  overpotential  as  function  of  total  current  density,  (b)  the 
fractional  coverage  of  catalytic  sites  by  hydrogen  and  (c)  the  fractional 
coverage  by  CO.  Notice  that  in  the  case  of  the  highest  CO  content  the 
curves  for  the  two  different  krs  values  almost  coincide. 

Without  the  RWGS  reaction  the  coverage  is  smaller,  simi¬ 
lar  to  what  was  found  in  the  case  without  CO  oxidation.  As 
the  current  density  and  the  overpotential  increase,  the  cov¬ 
erage  by  CO  in  both  cases  tends  to  the  same  value,  e.g.  the 
additional  CO2  poisoning  effect  disappears. 

3.3.  Effect  of  COi  in  H2,  finite  utilisation 

The  CO2  sensitivity  of  an  anode  can  be  tested  by  measur¬ 
ing  the  performance  of  a  cell  operated  in  FG  with  increasing 
amounts  of  CO2.  Here,  we  present  model  results  obtained 
for  a  gas  mixture  saturated  with  water  at  80  °C,  1.5  bar.  A 
H2  utilisation  of  67%  was  assumed.  To  enable  compari¬ 
son  with  recent  experimental  data  in  our  laboratory,  cell 
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Fig.  3.  Results  from  the  kinetic  model  obtained  with  pFF  =  0.8  bar,  pC02 
=  0.2  bar,  and  pCO  =  20  (jl  bar.  Electrochemical  oxidation  of  CO  was 
included  with  ktc  =  3  x  1 0  6  A/crrr .  The  dotted  lines  with  open  symbols 
correspond  to  the  case  kIS  =  0,  the  drawn  curves  to  kIS  =0.02  A/cm2 . 

performances  at  a  fixed  current  density  were  calculated. 
The  anode  was  modelled  with  the  parameter  values  from 
Table  1,  the  cathode  and  ohmic  losses  were  calculated  using 
the  parameter  values  from  Table  2.  The  values  were  obtained 
by  fitting  Eq.  (13)  to  an  experimental  V-I  characteristic  of 
a  cell  operated  under  the  above  conditions  with  Eb/air. 

In  Fig.  4  the  effect  of  increasing  CO2  content  is  presented. 
As  expected  the  performance  decreases  with  increasing  CO2 
fraction.  However,  as  Fig.  4b  shows,  the  amount  of  CO  in 
the  anode  outlet  gas  decreases  with  increasing  CO2  content. 
The  main  reason  for  this  effect  is  that  the  total  flow  increases 
when  the  CO2  fraction  increases,  as  a  result  of  the  fixed 
hydrogen  utilisation.  When  there  is  no  oxidation  of  CO  the 
total  amount  formed  in  the  RWGS  reaction  can  easily  be 
calculated  from  the  outlet  concentration.  Fig.  4c  shows  that 
this  amount  increases  with  the  CO2  fraction,  as  expected. 
The  results  can  be  compared  to  those  obtained  with  a  mix¬ 
ture  of  H2/N2/CO2  where  the  H2  flow  is  constant  but  N2  was 
gradually  replaced  by  CO2.  In  that  case  the  flow  is  constant, 
but  larger  than  the  H2/CO2  flow  due  to  the  extra  nitrogen. 
The  CO  concentration  at  the  anode  outlet  increases  slightly 


Table  2 

Values  for  the  parameters  in  Vst  (ohmic  and  cathode  losses). 


Ki  (V) 

0.756 

K2  (V) 

0.024 

(ohm  cm2) 

0.173 

Ka  (V) 

9.29  x  10~4 

K5  (cm2/A) 

3.1888 

Fig.  4.  Cell  performance,  CO  concentration  in  the  anode  outlet  and  rate 
of  CO  formation  as  a  function  of  the  CO2  content  in  FIs  (or  H2/N2).  The 
values  were  obtained  at  a  current  density  of  350  mA/cm2 ,  other  conditions 
see  text.  The  base  case  (♦)  was  calculated  using  the  parameter  values  of 
Table  1  with  kcc  =  0  A/cm2  and  krs  =  0.02  A/cm2  bar  and  using  FI2/CO2 
as  feed.  The  dashed  curves  with  open  symbols  (Oj  were  calculated  with 
50%  FF,  the  indicated  fraction  of  CO2,  balance  N2.  The  other  curves 
were  calculated  as  the  base  case  but  with,  respectively,  100  x  kIS  (A), 
0.1  x  kd  (*),  2  x  kic/kdc  (O)  and  ltc  =  3  x  10-6A/cm2  (I). 


as  result  of  the  larger  pCOi  but  is  smaller  than  in  the  case 
of  H2/CO2  when  the  flow  is  smaller.  The  actual  amount  of 
CO  formed  is  larger  than  in  H2/CO2,  in  spite  of  the  smaller 
fraction  of  hydrogen.  This  means  that  the  larger  CO  con¬ 
centration  in  H2/CO2  inhibits  further  reaction  of  CO2. 

With  respect  to  Fig.  4b  it  must  be  noted  that  the  concen¬ 
trations  shown  here  can  be  larger  than  what  is  expected  for 
thermodynamic  equilibrium.  In  fact  the  fuel  cell  conditions 
are  not  equilibrium  conditions,  and  the  rate  constants  of  the 
reactions  have  been  chosen  such  that  they  reflect  certain  as¬ 
sumptions  concerning  the  kinetics.  The  values  chosen  here 
imply  that  the  normal  water  gas  shift  reaction  cannot  take 
place,  the  activation  of  water  needs  a  certain  overpotential. 
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When  the  rate  constant  of  the  RWGS  reaction  is  increased 
by  a  factor  of  100,  the  performance  drops  (Fig.  4a).  However, 
Fig.  4b  shows  that  the  amount  of  CO  at  the  outlet  is  not  much 
affected.  This  means  that  similar  amounts  of  CO  are  being 
formed  as  is  confirmed  by  Fig.  4c.  As  was  already  shown 
above,  a  higher  rate  constant  results  in  a  larger  value  for 
Oco  in  the  stationary  state,  and  this  results  in  a  lower  value 
for  0\\.  This  again  reduces  the  amount  of  CO  formed.  Since 
the  actual  rate  of  the  RWGS  reaction  (given  by  Eq.  (14)) 
depends  through  On  in  a  complex  way  on  the  rate  constant 
£rs,  the  amount  of  CO  in  the  outlet  is  not  a  direct  measure 
for  the  rate  constant  of  the  of  the  RWGS  reaction. 

Another  parameter  that  will  have  a  large  effect  on  the 
CO  coverage  is  the  ratio  kiJkac.  Increase  by  a  factor  of  two 
results  in  a  decrease  of  the  maximum  value  for  Oco  at,  e.g. 
10  ppm  from  0.77  to  0.62.  The  effect  is  an  increase  of  the  cell 
performance  (Fig.  4),  as  lower  Oco  results  in  higher  values 
for  On.  This  also  increases  the  rate  of  CO  formation  in  the 
RWGS.  The  CO  content  at  the  outlet  is  therefore  larger  than 
in  the  case  of  small  k&Jkac  but  the  net  effect  is  an  increase 
of  performance. 

Also  the  dissociation  rate  of  hydrogen  will  influence  the 
CO2  poisoning  effect.  A  smaller  value  for  ka  results  in  an 
increased  performance  loss.  The  main  effect  of  a  lower  ka  is 
that  On  is  reduced.  This  results  in  a  lower  rate  of  formation  of 
CO  and  lower  CO  content  at  the  outlet,  but  not  in  a  reduction 
of  Oco-  Notice,  that  the  effect  of  this  lower  ka  on  the  cell 
performance  is  negligible  when  there  is  no  CO2  present.  It 
is  smaller  than  the  effect  of  replacing  50%  of  the  H2  by  N2 
as  Fig.  4a  shows. 

The  effect  of  CO  poisoning  can  be  mitigated  if  electro¬ 
oxidation  of  CO  is  possible,  as  shown  in  Fig.  4a.  The  amount 
of  CO  at  the  outlet  also  decreases  as  does  the  amount  of  CO 
formed  minus  oxidised  with  respect  to  the  base  case. 

3.4.  Reformate  feed,  finite  utilisation 

In  order  to  analyse  losses  with  a  typical  reformate  feed 
stream  containing  40%  H2,  40%  N2,  20%  CO2  and  10 ppm 
CO,  water  saturated  at  80  °C,  and  with  a  pressure  of  1.5  bar, 
several  V—I  curves  were  calculated  with  different  anode 
feeds.  The  ohmic  and  cathode  losses  were  calculated  accord¬ 
ing  to  Eq.  (13),  using  the  parameter  values  listed  in  Table  2. 
The  anode  losses  were  calculated  with  the  kinetic  model  and 
parameter  values  from  Table  1  (kec  =  3  x  10-6  A/cm2,  kIS 
=  0.02  A/cm2).  An  overall  utilisation  of  90%  was  assumed. 
Fig.  5  shows  the  results.  The  difference  between  the  curves 
obtained  with  pure  H2  and  with  H2/N2  represents  the  losses 
due  to  dilution  and  utilisation. 

When  part  of  the  nitrogen  is  replaced  by  CO2,  the  cell 
performance  decreases  due  to  catalyst  surface  blocking  as 
a  result  of  the  RWGS  reaction  and  CO  formation.  With  the 
parameter  values  chosen  here  CO2  poisoning  presents  an  ad¬ 
ditional  loss  comparable  to  the  dilution  and  utilisation  loss. 
Moreover,  at  current  densities  below  0.5  A/cm2  the  CO2  poi¬ 
soning  is  comparable  to  the  losses  obtained  in  a  feed  with- 


Current  Density  (A/cm2) 

Fig.  5.  V-I  curves  calculated  with  different  fuel  compositions.  Utilisation 
in  all  cases  90%.  For  details  of  the  calculation  see  text.  The  numbers  of 
in  legend  refer  to  the  fraction  of  the  gas  species  in  the  total  (dry)  flow. 

out  CO2  but  with  10  ppm  CO.  At.  At  higher  current  density 
the  rate  of  CO2  reduction  is  reduced  due  to  lower  amount 
of  H-adsorbed  sites  available  for  the  CO2  reduction  reac¬ 
tion,  hence  at  higher  current  density  the  effects  of  10  ppm 
CO  become  larger  than  those  of  CO2. 

In  Fig.  6,  the  CO  content  in  both  cases  is  shown  calcu¬ 
lated  at  a  cell  voltage  of  0.6  V.  In  the  case,  of  no  CO  in  the 
feed  there  is  a  gradual  increase  of  the  CO  content  in  the  an¬ 
ode  channel,  until  it  decreases  near  the  channel  exit  where 
the  overpotential  will  be  highest  and  CO  electro-oxidation 
will  be  faster  than  CO2  reduction.  In  the  case  of  no  CO2 
in  the  feed  the  CO  content  drops  gradually  from  the  initial 
value  as  a  result  from  electro-oxidation.  Notice  that  without 
electro-oxidation  the  CO  content  would  increase  also  in  this 
case. 

When  both  CO2  and  CO  are  present  in  the  feed,  the  CO 
content  is  on  the  whole  larger  than  in  the  cases  without  CO  or 
without  CO2  and  goes  through  a  maximum.  In  the  V—I  curves 
at  intermediate  current  density  the  effect  of  both  CO  and  CO2 
in  the  feed  is  an  additional  performance  loss,  but  the  effect 
is  not  additive.  Moreover,  at  high  current  densities  the  CO 
effect  seems  to  prevail,  whereas  at  low  current  densities  the 
CO2  effect  dominates.  This  is  similar  to  results  obtained  with 


Position  in  Channel 

Fig.  6.  The  CO  content  as  a  function  of  the  position  in  the  gas  channel 
in  the  case  of  20%  CO2  in  the  dry  feed  but  no  CO  (O),  no  CO2  but 
10  ppm  CO  (□),  10  ppm  CO  and  20%  CO2  in  the  feed  (A).  The  data 
correspond  to  point  Vceii  =  0.6  V  in  the  previous  figure. 
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the  kinetic  model  and  displayed  in  Fig.  3.  In  the  presence  of 
CO  the  H-coverage  of  the  catalyst  surface  area  is  small  and 
the  RWGS  reaction  is  suppressed. 

4.  Discussion 

The  above  analyses  based  on  a  simple  model  have  shown 
how  CO2  poisoning  effects  depend  on  the  properties  of  the 
anode  catalyst  material.  The  most  important  conclusions 
are  that  increase  of  the  rate  constant  of  the  RWGS  reaction 
and  reduction  of  the  hydrogen  dissociation  rate  enhance 
CO2  poisoning  effects,  whereas  CO  desorption  and  CO 
electro-oxidation  mitigate  the  effect.  The  amount  of  CO  at 
the  anode  outlet  does  not  correlate  with  the  degree  of  CO2 
poisoning.  The  main  cause  of  CO2  poisoning  is  that  there 
is  a  stationary  high  coverage  of  the  catalyst  surface  area  by 
reduced  CO2,  e.g.  CO.  This  reduces  the  amount  of  surface 
area  covered  by  hydrogen,  which  reduces  both  the  hydro¬ 
gen  oxidation  rate  as  well  as  the  CO  formation  rate  in  the 
RWGS  reaction. 

Work  in  our  own  laboratory,  as  well  as  work  by  others 
[11,12]  has  shown  that  different  catalyst  materials  do  indeed 
respond  differently  to  CO2.  Some  results  recently  obtained 
in  our  laboratory  are  shown  in  Fig.  7.  Notice  that  the  shape  of 
the  measured  curves  resembles  well  those  of  the  ones  calcu¬ 
lated  using  the  model.  Moreover,  the  observed  anode  losses 
are  in  the  range  0-100  mV,  comparing  well  to  those  obtained 
with  the  range  of  values  used  for  the  model  calculations  (see 
Fig.  4).  However,  given  the  simplicity  of  the  model  and  the 
number  of  parameters  affecting  the  effect  of  CO2  poisoning 
no  attempt  will  be  made  here  to  fit  the  parameter  values  of 
the  model  to  the  experimental  curves,  and  only  a  qualitative 
interpretation  of  the  observed  curves  will  be  given. 

Fig.  7  shows  that  PtRu  alloy  is  more  CO2  tolerant  than  Pt. 
It  is  well  established  that  PtRu  alloy  is  also  more  CO  tolerant 
than  Pt.  At  low  overpotential  (<0.2  eV)  this  is  not  ascribed 
to  enhanced  CO  oxidation  but  to  reduced  adsorption  due  to 
the  ligand  effect  [15,16].  According  to  the  model  analyses 
presented  here  this  would  explain  a  higher  CO2  tolerance.  In 
this  case  the  expected  CO  content  in  the  anode  outlet  would 


Fig.  7.  CO2  poisoning  effect  for  various  anodes  using  a  H2/CO2  feed.  The 
figure  shows  the  anode  losses  at  a  fixed  current  density  of  35()mA/cm2 
as  a  function  of  CO2  content.  Hydrogen  stoichiometry  1.5,  80  °C,  1.5  bar 
gas  pressure. 


be  larger  than  on  Pt.  However,  it  is  also  quite  possible  that 
the  rate  constant  of  the  RWGS  reaction  is  lower  on  PtRu. 
Less  interaction  with  CO  may  also  imply  less  interaction 
with  CO2  and  intermediates.  The  final  source  of  increased 
CO2  tolerance  would  be  an  increased  hydrogen  dissociation 
rate.  However,  Gasteiger  et  al.  [17]  found  that  on  Pt  and 
PtRu  the  electro-oxidation  reaction  of  H2  was  equally  fast, 
and  much  slower  on  pure  Ru,  so  it  seems  this  option  must 
be  discarded. 

PtMo  catalysts  have  a  good  CO  tolerance,  which  is  as¬ 
cribed  to  CO  electro-oxidation  at  low  (<0.2  V)  overpoten¬ 
tials  [18].  On  the  other  hand,  the  CO2  tolerance  is  not  very 
good.  Ball  et  al.  [11]  reported  that  the  CO2  tolerance  de¬ 
creased  in  the  order  PtRu,  Pt,  PtMo,  as  corroborated  by  our 
results  (Fig.  7).  Recently,  Urian  et  al.  [12]  presented  more 
extensive  data  on  CO2  poisoning  of  PtRu,  Pt,  and  PtMo  cat¬ 
alysts,  which  showed  that  at  current  densities  larger  than 
0.5  A/cm2  Pt  becomes  the  most  CO2  intolerant  species  of 
the  three.  The  structure  of  the  PtMo  catalyst  surface  will  de¬ 
pend  strongly  on  the  composition  as  well  as  the  preparation 
method,  which  explains  differences  in  CO  and  CO2  toler¬ 
ance  observed  by  different  workers.  However,  in  general  it 
is  assumed  that  Mo  particles  at  the  surface  of  both  PtMo  al¬ 
loys  and  bimetallics  can  easily  form  (hydro)oxides,  which 
accounts  for  the  easy  activation  of  water,  a  prerequisite  for 
the  electro-oxidation  of  CO  at  low  overpotential.  But  the 
presence  of  oxidised  Mo  blocks  Pt  sites  for  hydrogen  dis¬ 
sociation  [18].  Translated  to  the  model  presented  here  this 
leads  to  reduced  values  for  the  hydrogen  dissociation  rate 
kd,  which  has  a  negative  effect  on  the  CO2  tolerance  as  illus¬ 
trated  by  Fig.  4.  The  observation  made  by  Urian  et  al.  [12] 
that  at  high  current  densities  the  CO2  tolerance  increases 
with  decreasing  Pt:Mo  ratio  may  be  related  to  the  oxidation 
of  CO2  reduction  species  facilitated  by  the  presence  of  Mo. 
The  kinetics  of  CO2  reduction  on  supported  Pt  and  PtMo 
catalysts  are  at  present  being  studied  in  our  laboratory  and 
will  be  the  subject  of  a  subsequent  publication. 

In  the  literature  it  has  been  observed  that  the  CO2  tol¬ 
erance  of  anodes  also  seems  to  depend  on  the  microstruc¬ 
ture.  It  was  concluded  that  a  large  fraction  of  electrochemi- 
cally  inactive  particles,  i.e.  particles  not  covered  by  Nafion, 
increased  the  CO2  intolerance  [5,19].  An  explanation  sug¬ 
gested  was  that  electrochemically  inactive  particles  are  more 
active  in  CO  formation  as  the  amount  of  adsorbed  hydrogen 
would  be  larger  on  these  particles.  This  CO  would  than  des¬ 
orb  into  the  gas  phase  and  poison  active  particles.  In  ref.  [3] 
it  was  shown  that  electrochemically  inactive  particles  are  ac¬ 
tive  catalysts  for  CO  formation  from  H2  and  CO2.  However, 
as  Fig.  1  shows,  at  least  in  the  case  of  a  high  kd  the  values 
of  %  and  7(  0  do  not  depend  much  on  the  current  density, 
i.e.  on  whether  particles  are  electrochemically  active  or  not. 
This  means  that  an  increased  effect  of  CO  formation  and 
desorption  into  the  gas  phase  from  inactive  particles  is  not 
expected.  The  analyses  presented  here  rather  suggest  that  the 
reduced  values  of  the  hydrogen  dissociation  rate  ka,  which 
result  for  an  electrode  with  an  increased  fraction  of  inactive 
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particles,  is  the  reason  for  the  decrease  in  CO2  tolerance.  The 
data  presented  by  Ralph  and  Hogarth  [19]  show  that  the  ac¬ 
tive  surface  area  in  various  anodes  with  similar  loading  can 
vary  by  a  factor  of  six,  which  would  correspond  to  similar 
variations  of  kd.  As  Fig.  4  shows,  a  reduction  of  ka  results 
in  reduced  CO2  tolerance  as  well  as  reduced  CO  formation. 

The  results  have  also  shown  that  in  true  reformate,  where 
CO  is  present  as  well  as  CO2,  that  CO2  will  present  addi¬ 
tional  losses,  but  the  losses  due  to  CO  and  CO2  are  not  ad¬ 
ditive.  At  high  CO  content  CO2  reduction  cannot  take  place 
as  there  is  not  sufficient  adsorbed  hydrogen  present  on  the 
catalyst. 

5.  Conclusions 

Model  calculations  have  been  carried  out  to  study  the 
CO2  poisoning  effect  on  anode  catalyst  materials  of  PEM 
fuel  cells.  The  main  effect  of  CO2  poisoning  is  that  due  to 
the  RWGS  reaction  a  large  part  of  the  catalytic  surface  area 
becomes  inactive  for  H2  dissociation.  Subsequent  desorp¬ 
tion  of  CO  from  the  catalyst  surface,  transport  down  the  gas 
channel,  and  subsequent  re-adsorption  of  CO  plays  a  minor 
role.  The  main  reason  for  this  is  that  a  large  blockage  of  the 
surface  area  inhibits  further  formation  of  CO  in  the  RWGS 
reaction. 

It  was  found  that  a  high  rate  constant  of  the  RWGS  re¬ 
action  increases  the  anode  polarisation  losses,  as  does  a  re¬ 
duced  rate  constant  of  the  hydrogen  dissociation  reaction. 
The  effects  are  mitigated  by  a  high  ratio  of  the  CO  desorp¬ 
tion  and  adsorption  rate  constants,  as  well  as  by  a  high  CO 
electro-oxidation  rate  constant. 

In  true  reformate  where  CO  is  present  besides  CO2,  CO2 
poisoning  will  have  the  largest  effect  when  the  CO  content 
is  small  and  at  relatively  low  current  density. 
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